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Abstract 

We describe the geometry of geodesies on a Lorentz ellipsoid: give 
explicit formulas for the first integrals (pseudo-confocal coordinates), 
curvature, geodesically equivalent Riemannian metric, the invariant 
area-forms on the time- and space-like geodesies and invariant 1-form 
on the space of null geodesies. We prove a Poncelet-type theorem 
for null geodesies on the ellipsoid: if such a geodesic close up after 
several oscillations in the "pseudo-Riemannian belt", so do all other 
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g : 1 Introduction 
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The geodesic flow on the ellipsoid in Euclidean space is a classical exam- 
ple of a completely integrable dynamical system whose study goes back to 
Jacobi and Chasles. We refer to [3] [13], [15] for a modern treatment of this 
subject; see also [U E] [9] [10], [TT] [22] for a sampler of recent work. In a 
recent paper [12] . we considered ellipsoids in pseudo- Euclidean spaces of ar- 
bitrary signatures and extended, with appropriate adjustments, the theorem 
on complete integrability of the geodesic flow to this case. We also defined 
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pseudo-Euclidean billiards and proved that the billiard map inside an ellip- 
soid in pseudo-Euclidean space is completely integrable. 

This note is devoted to the "case study" of geodesies on an ellipsoid in 
three dimensional space with the metric dx 2 + dy 2 — dz 2 . Recall that the 
pseudo-sphere x 2 + y 2 — z 2 = —1 in Minkowski space provides a famous 
example of a Riemannian metric of constant negative curvature, a model 
of the hyperbolic plane, and this gives a further motivation to investigate 
quadratic surfaces in Minkowski space. 

We shall study an ellipsoid 



and we assume that the general position condition a > b holds. The induced 
metric ( , ) on the ellipsoid degenerates along the two curves 



that will be referred to as the "tropics" . The induced metric is Riemannian 
in the "polar caps" and Lorentz in the "equatorial belt" bounded by the 
tropics. We shall see in Section [6] that the Gauss curvature of the ellipsoid 
is everywhere negative (it equals — oo on the tropics). At every point of the 
equatorial belt, one has two null directions of the Lorentz metric; on the 
tropics, these directions merge together. 

Every geodesic curve j(t) on the ellipsoid ([1]) is of one of the three types: 
space-like (positive energy (7', 7') > 0), time-like (negative energy (7', 77 < 
0) or light-like (zero energy or null (7', 7') = 0). The last two types exist 
only in the equatorial belt. In fact, a whole new phenomenon which we are 
dealing with in the Lorentz compared to the Euclidean one, is the 

presence of null geodesies, which "separate" the space-like and time-like ones. 
This makes the pseudo-Riemannian geometry of the problem quite peculiar, 
and rather different from its Riemannian counterpart. 

Let us summarize some relevant results from [T2], specialized to the 
three dimensional case. Include the ellipsoid into a pseudo-confocal family of 
quadrics M\ given by the equation 




a, 6, c > 0, 



(1) 




(2) 



x 2 y 2 z 2 



1. 



(3) 



a + A b + A c — A 



It is shown in [12] that: 
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1. Through every generic point Q(x,y,z) in space, there pass either three 
or one pseudo-confocal quadrics. In the former case, two of the quadrics 
have the same topological type, and the quadrics are pairwise orthog- 
onal at point Q. If Q is a generic point of the ellipsoid ([I]) then there 
exists another pseudo-confocal ellipsoid and a pseudo-confocal hyper- 
boloid of one sheet, passing through x. 

2. A generic space- or time- like line I is tangent to either two or no pseudo- 
confocal quadrics. In the former case, the tangent planes to these 
quadrics at the tangency points with I are pairwise orthogonal. In 
particular, a generic line, tangent to the ellipsoid (jTJ), is tangent to 
another pseudo-confocal quadric. 

3. The tangent lines to a given space-like or time-like geodesic on the 
ellipsoid ([1]) forever remain tangent to a fixed pseudo-confocal quadric; 
the respective value of A in ([3]) can be considered as an integral of the 
geodesic flow. 

The spaces of space- or time-like lines in pseudo-Euclidean space, and 
more generally, the spaces of space- or time-like geodesies in a pseudo- 
Riemannian manifold, carry symplectic structures. These structures are ob- 
tained by symplectic reduction: restrict the canonical symplectic structure 
of the (co)tangent bundle to the unit energy hypersurface (v,v) = ±1 and 
quotient out the one-dimensional kernel of this restriction. In the case of a 
pseudo-Riemannian (Lorentz) surface, we obtain an area form on the set of 
space- and time-like geodesicsQ 

Now, to billiards. In general, the billiard dynamical system in a pseudo- 
Riemannian manifold with a smooth boundary describes the motion of a free 
mass-point ("billiard ball"). The point moves along a geodesic with constant 
energy until it hits the boundary where the elastic reflection occurs: the 
normal component of the velocity instantaneously changes sign whereas the 
tangential component remains the same, see [2U1 EI] for general information. 
The billiard reflection is not defined at the points where the normal vector 
is tangent to the boundary. The billiard ball map acts on oriented geodesies 
and takes the incoming trajectory of the billiard ball to the outgoing one. 

1 It is shown in [T^] that the space of light-like geodesies has a contact structure; it is 
trivial in the case of a surface, as light-like geodesies form a one-dimensional set, and we 
shall not use this structure in the present paper. 
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This map preserves the type of a geodesic (space-, time-, or light-like). Fur- 
thermore, the billiard ball map, acting on the space- or time-like geodesies, 
preserves the symplectic structure (area form, in two-dimensional case) de- 
scribed aboven 




It is proved in [12] that the billiard ball map with respect to a quadric 
in Minkowski space is integrable as well: if an incoming space- or time- 
like billiard trajectory is tangent to two pseudo-confocal quadrics then the 
outgoing trajectory remains tangent to the same two quadrics. 

In general, complete integrability of a billiard system implies configura- 
tion theorems for closed billiard trajectories, see, e.g., [20J, [2Tj. In particular, 
the integrability of the billiard inside an ellipse implies the classical Poncelet 
Porism depicted in figure [TJ let 7 C T be two nested ellipses and let a point 
of T be a vertex of an n-gon inscribed in T and circumscribed about 7; then 
every point of T is a vertex of such an n-gon, see [8]. 



In Section [5] we prove a Poncelet- style closure theorem involving null 
geodesies on the ellipsoid: if some chain of such geodesies in the equatorial 
belt closes up after n oscillations between the tropics then so does every such 
chain. 

It is proved in [T2| that the billiard inside an ellipse in the Lorentz plane 
is completely integrable as well. Restricted to light-like lines, the billiard 
ball map in an oval becomes a circle map depending only on the two null 
directions. In Section [7| we discuss the problem when this map is conjugated 
to a rotation. 

2 And, acting on light-like geodesies, preserves the contact structure. 





Figure 1: Poncelet Porism 
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2 Joachimsthal integral 



The geodesic flow on the triaxial ellipsoid in Euclidean space possesses the 
classical Joachimsthal integral. In this section we describe its pseudo- Euclidean 
version. This integral says analytically that the lines tangent to a geodesic 
curve on the ellipsoid are tangent to a fixed pseudo-confocal quadric: in this 
dimension, all first integrals of the geodesic flow, restricted to a constant 
energy hypersurface, are functionally dependent. 

Choose a Minkowski normal at point (x, y, z) of the ellipsoid ([I]) as follows: 

\rf ^ fx y z 
N{x,y,z) = -,-r, — 
\a b c 

Then 

x 2 y 2 z 2 

(N(x, y, z), N(x, y,z)) = — + — - —. 

a o c 

We see that the normal is space-like in the equatorial belt, time-like in the 
polar caps and null on the tropics. 

Let (u, v, w) denote a tangent vector to the ellipsoid. 

Proposition 2.1 The following function is an integral of the geodesic flow: 

x 2 y 2 z 2 \ ( u 2 v 2 w 2 



Proof. Parameterized geodesies j(t) are extrema of the functional 

( 7 (t), 7 (t)) dt- 



they satisfy the Euler-Lagrange equation j(t) = A(t)iV(7(i)) where X(t) is a 
Lagrange multiplier. 

Let 7 (t) be a geodesic with 7(t) = (u, v,w). Then 7 (t) = \(t)N(y(t)), 
and hence 

Am m),N(y(t))) 

U <JV( T (t)),JV( T (t))>- 
Since ( 7 (t), N(y(t))) = 0, we have: 

,2 „,2 



m),N(y(t))) = -(i(t),N(l(t))} = -((u,v,w), (-,p--)> = - ( - + 

\a o c J \ a o 



and it follows that 

u 2 I v 2 I ui 2 

A « = ~ kA_k < 5 » 

a 2 "T b 2 c 2 

(the denominator vanishes precisely on the tropics). 
We want to prove that J = 0. Indeed, 

fxu yv zw\ ( u 2 v 2 w 2 \ ( x 2 y 2 z 2 \ f uu vv ww 



\ a 2 b 2 c 2 ) \ a b ' c J ' \a 2 ' b 2 c 2 J \ a h <- 

rxu ^yv zw\ fu 2 v 2 w2 \_^y^ ( x2 V 2 2 \ ( xu _|_ V v zw 
\a 2 b 2 c 2 ) V a b c J \a 2 b 2 c 2 J V a 2 b 2 c 2 



the last equality due to (JSJ). □ 



3 Describing the geodesies 
3.1 Global behavior of geodesies 

As we mentioned earlier, the metric in the polar caps is Riemannian, while 
in the equatorial belt it is Lorentzian. A geodesic in a polar cap crosses it 
from tropic to tropic. 

Consider the equatorial belt. The light-like geodesies traverse the belt 
from one tropic to another. There are two kinds of them depending on 
the slope, positive or negative; we call these null geodesies right and left, 
respectively. The time-like geodesies are squeezed between like-like ones and 
hence also go from tropic to tropic. 

Let us discuss the behavior of space-like geodesies. The equator is a closed 
geodesic. Figure [2] depicts two other space-like geodesies (first two pictures); 
the first one intersects the equator and proceeds to the other tropic, and 
the other geodesic starts at a tropic and returns to the same tropic without 
reaching the equator. These two geodesies represent generic behavior as the 
next lemma shows. 

Lemma 3.1 Let j(t) = (x(t),y(t), z(t)) be a geodesic in the Northern part 
of the equatorial belt, i.e., z{t) > 0. Then z{t) > 0. In the Southern part of 
the equatorial belt, z(t) < 
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Figure 2: Space-like geodesies in the equatorial belt 



Proof. Since j(t) = \(t)N(y(t)), we have: z(t) = —\{t)z{t) / c. According 
to (jSJ), A < in the equatorial belt, hence z(t) > for z > (i.e., in the 
Northern part) and z(t) < for z < (in the Southern part). □ 

It follows that a geodesic in the Northern part of the equatorial belt is 
convex downwards, and in the Southern part of the equatorial belt is convex 
upwards. It also follows from the proof that the equator, as a closed geodesic, 
is exponentially unstable. 

An intermediate position between the two types of generic space-like 
geodesies is occupied by the geodesies that start at, say, Northern tropic 
and monotonically descend to the equator having the latter as a limit cycle, 
see figure IU on the right. 

3.2 Local behavior of geodesies near the tropics 

The behavior of geodesies near tropics is described by the following 

Proposition 3.2 A geodesic on the ellipsoid (QJ) may reach a tropic only in 
the null direction. 
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Proof. For light-like geodesies, there is nothing to prove, and a time-like 
geodesic is confined in the wedge between two null directions which merge 
together on a tropic. It remains to consider unit energy space-like geodesies. 
According to Proposition 12.11 one has, along a geodesic: 



As the point approaches a tropic, the first factor on the left hand side of 
the first equality goes to zero, and hence the second factor blows up. The 
direction of the geodesic does not change if we rescale the tangent vector: 



maintaining u 2 + v 2 + w 2 = 0(1). As (x, y, z) approaches the tropic, we have 

u 2 + v 2 - w 2 = /i 2 -y 0, 

therefore the geodesic has the null direction. □ 

It follows from Proposition 13.21 that the null lines tangent to the ellipsoid 
along the tropics are tangent to infinitely many pseudo-confocal quadrics. 
Indeed, infinitely many geodesic hit the tropics at any given point in the 
null direction, and the tangent lines to each geodesic are tangent to a fixed 
pseudo-confocal quadric, see Section [TJ In contrast, we have the following 
proposition which provides an alternative proof of Proposition 13.21 

Proposition IUT2T Let P be a point of a tropic and £ a space-like line tangent 
to the ellipsoid at P. Then £ is not tangent to any pseudo-confocal quadric. 

Proof. Assume that £ is tangent to M\ at point Q. Denote by N a normal 
vector to Mq at P and by r] a normal vector to M\ at Q. Then N lies in the 
tangent plane TpM . The restriction of the ambient metric to this plane is 
degenerate, and N spans the kernel of this restriction. 

According to the general theory, see Section [H rj is orthogonal to N, hence 
rj lies in N 1 - = TpM^. Since rj is also orthogonal to i, the normals r) and iV 
are collinear. Therefore the tangent planes TpMo and TqM\ coincide. 

The projection of pseudo-confocal family along the line £ is a pseudo- 
confocal family of conies in the Lorentz plane, see [12]. It follows that the 




) 



const, u 2 + v 2 — w 2 = 1. 



U, V, w 
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Figure 3: Pseudo-confocal family of conies 4- _2L_ = 1 

projections of M Q and M\ are tangent to each other. However no two conies 
in a pseudo-confocal family are tangent (see figure [3]), a contradiction. □ 

Proposition 13.21 is a manifestation of a more general phenomenon which 
holds for Lorentz surfaces. Suppose, on an analytic Lorentz surface S, the 
Lorentz metric changes into Riemannian (i.e., the two null direction coincide) 
along a smooth curve T. We call a point of V regular if the null direction at 
this point is transversal to the curve T itself. 

Proposition 3.3 In a neighborhood of a regular point the Lorentz metric on 
the surface S is conformally equivalent to the metric dy 2 —xdx 2 in a neighbor- 
hood of the origin. The null geodesies on S near this point are diffeomorphic 
to the family of cusps y = x 3//2 + C , see figure^ 

Proof. Let a(x, y)dx 2 + b(x, y)dxdy + c(x, y)dy 2 be a Lorentz metric on S. 
Then the equation of null geodesies on 5 is a(x, y)+b(x, y)y'+c(x, y){y') 2 = 0, 
an implicit differential equation of the first order. The condition of transver- 
sality at a regular point of T means that we consider a so-called "regular 
singular point" of the implicit differential equation F(x,y,y') = 0. Accord- 
ing to the theorem of Cibrario (see a discussion in [2]) the normal form of 
this differential equation at a regular singular point is (y') 2 = x, which is the 
equation of null geodesies in the metric dy 2 — xdx 2 . Thus the null geodesies 
is the family of cusps described above. 
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Figure 4: Family of cusps 



The family of null geodesies fixes the conformal class of the Lorentz met- 
ric. Due to analyticity, this metric uniquely extends beyond the curve T, into 
the Riemannian domain of the neighborhood of the regular point. □ 

Remark 3.4 Note that every geodesic for the metric dy 2 — xdx 2 which hits 
the "tropic" T = {x = 0}, does it with a horizontal velocity. Indeed, consider 
a non-vertical space-like geodesic with the velocity vector (u, v) of the unit 
Lorentz length: v 2 — xu 2 = 1. Note that v 2 ^ 1 for any x ^ 0. Since the 
metric is invariant with respect to y-translations, the velocity vector along 
any geodesic conserves its ^-component (in addition to the conservation of 
its Lorentz length). Then, at the moment of "impact" with the tropic, the 
■u-component has to become infinite: u 2 = (v 2 — l)/x — > oo as x — > 0, that 
is, the velocity becomes horizontal on the tropic. 

Although, for a metric which is only conformally equivalent to this normal 
form, one cannot use the invariance of the v -component, one can bound it 
above and below, and hence the u-component has to become infinite anyway, 
i.e., the same conclusion on the horizontality of the velocity holds, due to its 
robustness. This consideration implies Proposition 13.21 for the ellipsoid, once 
one checks that the null directions are everywhere transversal to the tropics, 
which is straightforward. 
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3.3 Intersections with pseudo-confocal quadrics 

The topology of a pseudo-confocal quadric (j3J) depends on the position of A 
relative the three numbers —a, —b and c: if A < —a then M\ is a hyperboloid 
of two sheets, if —a < A < — b then M\ is a hyperboloid of one sheet, if 
— b < A < c then Mx is an ellipsoid, and if c < A then M\ is again a 
hyperboloid of one sheet. Only the second and the third kinds intersect the 
original ellipsoid, M . See figure [5] for all five quadrics and figure [6] for an 
ellipsoid M\ with — b < A < c intersecting M . 




Figure 5: Four pseudo-confocal quadrics, along with the initial ellipsoid 
The intersection curve T\ = M fl M x is given by the equation 

" 2 +ItAt = ^At- (6) 



a(a + X) b(b + X) c(c - A) 

In the limit A —>■ 0, this becomes the equation of the tropic (j2j). As A —>■ c, 
the curve T\ tends to the equator. The projection of the intersection curve 
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Figure 6: A pseudo-confocal ellipsoid intersecting the initial one 

T A on the (x, y)-plane is a conic 

a + c 9 b + c 9 
r -I u z — 1 

a(a + X) b(b + X) y 

If —a < A < — b then this conic is a hyperbola, and if — b < A < c it is an 
ellipse. 

Fix a value of A. Given a generic point P of the ellipsoid Mo, the number 
of tangent lines from P to the conic TpMo H Ma may equal 2 or 0, and 
the curve T\ separates these two domains, say, U\ and V\ (the number of 
tangents equals 1 for P G 

Consider a geodesic 7 on M . The lines tangent to 7 are tangent to 
some pseudo-confocal quadric M\. Therefore 7 is confined to U\ and, at any 
point P G £/a ; the geodesic 7 may have only two directions: the tangent 
directions from P to the conic TpM n Ma. The geodesic 7 cannot intersect 
the boundary Ta but can touch it. 

Note that the equator plays a special role: no line tangent to the equator 
is tangent to any pseudo-confocal quadric (except M ). Indeed, such a line 
lies in the horizontal plane, and the trace of the pseudo-confocal family ([3]) 
in the horizontal plane is the confocal family of conies. But confocal conies 
have no common tangents. 

3.4 Reflection of geodesies off the tropics 

Define "reflection" of a geodesic 7 off a tropic as the geodesic tangent to the 
same pseudo-confocal quadric as 7 or, equivalently, having the same value 
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of the Joachimsthal integral (normalizing the energy to ±1). According to 
Proposition 13.21 a geodesic and a reflected one have the null direction at the 
impact point and hence are tangent to each other. 

Proposition 3.5 1). Let 7 be a geodesic in a polar cap tangent to a curve 
T\ = Mq fl M\, and let 71 be the reflection 0/7 in the tropic. Then 71 is also 
tangent to Y\ (see figure\2ty. 

2). Let 7 be a space-like geodesic in the equatorial belt that does not intersect 
the equator, and let 71 be the reflection of 7 in the tropic. Then 7 is tangent 
to a curve Y\, separating 7 from the equator, and 71 is also tangent to Y\ 
and therefore disjoint from the equator, see figure 



Figure 7: Reflected geodesies in a polar cap 



Proof. Since the geodesic 7 is tangent to T A , the tangent lines to 7 are 
tangent to M\. By definition, the tangent lines to 71 are tangent to M\ as 
well, and hence 71 will touch T\, and the first claim follows. 

Likewise, the tangent lines to a space-like geodesic 7 in the equatorial 
belt are tangent to a pseudo-confocal geodesic M\. Since 7 does not reach 
the equator, the curve T x separates the tropic and the equator. Arguing as 
above, the reflected geodesic 71 is also tangent to T x and therefore disjoint 
from the equator. □ 
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Figure 8: Reflected geodesies in the equatorial belt (unfolded) 

3.5 Geodesies on degenerate ellipsoids 

It is useful to visualise the degenerations of the ellipsoid (0Q) into a two-sided 
flat surface, "a pancake." First consider the limit c — > 0, under which the 
ellipsoid gets squeezed to an ellipse with semi-axes a and b in the Euclidean 
plane {z = 0}. Then the (Lorentz) equatorial belt disappears, whereas each 
(Riemannian) polar cap becomes the interior of an ellipse. The geodesies 
in the polar caps become straight lines in the limit, and their reflection off 
the tropics becomes the billiard reflection in the ellipse; the Joachimsthal 
integral describes the confocal ellipse to which a billiard trajectory remains 
tangent. 

Compare this with the other degeneration, b — > 0. In the latter case the 
ellipsoid becomes (the interior of) an ellipse in the Lorentz plane {y — 0}. 
The polar caps (and tropics) become squeezed to the arcs of the ellipse z > 
and z < in this plane. The equatorial belt becomes the double cover of 
the interior of the ellipse. The geodesies also become straight lines in the 
Lorentz plane. The null geodesies have two prescribed slopes in this plane: 
x = ±z. We will see that the dynamics of the null geodesies in the ellipsoid 
is an interesting extension of the corresponding Lorentz billiard inside the 
ellipse, restricted to oriented null lines. 
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4 Area form on the space of geodesies 



In this section we compute the area form u on the space of time-like geodesies, 
described in Section [TJ We also define a 1-form on the space of light-like 
geodesies which will play the central role in the next section. 

Let us characterize a time-like or a light-like geodesic by its intersection 
with the equator of the ellipsoid ([I]). The equator is parameterized as Q(t) = 
(y/acost, v^sint, 0) for t G R/27rZ. Let a geodesic make (Euclidean) angle 
a with the equator. Then (t, a) are coordinates in the space of geodesies 
intersecting the equator. The null directions correspond to a = tc/4 and 
37r/4, and for time-like geodesies, 7r/4 < a < 37r/4. We set 

f(t) = a/ a sin 2 1 + 6cos 2 t, r 



V tan 2 a — 1 



Note that for a = 7r/2, the value of r is well-defined: r = 0. The following 
proposition describes the Joachimsthal integral J and the area form u in 
terms of the (t, a) coordinates. 

Proposition 4.1 One has: 



J{t,a) = J { R , tu = f(t)drAdt. (7) 

Proof. For the tangent vector to the equator, one has: 

Q'(t) = {-^/Esmt,Vbcost,0), (Q' (t) , Q' (t)) = f{t). 

The geodesic corresponding to (t, a) has a tangent vector 

(—y/asmt, Vbcost, f(t) tana). 

If the geodesic is time-like, we normalize the tangent vector so that its squared 
length is —1: 

/ \ / — \fa~Tsmt Vbrcost \ 

( u ^=[ m '-7^T' rtana J- (8) 

To obtain J, substitute (u, v,w) from (jHJ) and 

(x, y, z) = (yfa cost, Vb sin t, 0) (9) 
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to fll]); this yields the first formula ([7J). 

Identify the cotangent and tangent bundles via the metric. Then the 
canonical symplectic form on the cotangent bundle becomes 

uj = du A dx + dv A dy — dw A dz. 

Formulas and describe a section of the tangent bundle, and the pull 
back of the symplectic structure uj is given by the second formula ([7]). □ 



As a consequence of Proposition 14.11 we define a natural 1-form on the 
space of null geodesies (this space consists of two disjoint circles correspond- 
ing to the right and left null geodesies). The set of null geodesies is given, in 
(t, a) coordinates, by a = ir/4 or 37r/4, which corresponds to r = oo. Note 
that both uj and J blow up as one approaches the space of null geodesies, 
that is, as r — > oo. However, their ratio is well-defined. 

Lemma-Definition 4.2 The 1-form h(t)dt, given by the condition that, in 
the limit t — > oo, 

d(J 1/2 ) A h{t)dt = uj, (10) 
is uniquelly defined. Explicitly, this equation holds for h{t) given by 



,/\ fit) / a sin 2 1 + b cos 2 1 /11N 

nit) = const ■ — - = const ■ \ ^ — ill) 

v / c + f 2 {t) V c + asin 2 t + 6cos 2 t 

with an appropriate constant factor. 

Proof. Up to a constant multiplier, one has: 

J l/2 = r^VW) (l + O (1)) , 

dJ 1 ' 2 = ^c + f 2 (t) dr + g{t, r)dt + (J^J 
for some function g. Therefore 

as needed. □ 



hence 



16 



5 Poncelet-style closure theorem 



Define a map T of the equator to itself as follows. Given a point P of the 
equator, consider the right null geodesic through P until it intersects the 
Northern tropic at point Q, and then consider the left null geodesic through 
Q until it intersects the equator at point Pi. Set: T(P) = P\, see figure El 



Q 




Figure 9: The right null geodesic through P and the left null geodesic through 
Pi hit the Northern tropic at the same point Q. 

Theorem 5.1 Suppose that some point P of the equator is k-periodic, that 
is, T h (P) = P for some positive integer k. Then every point of the equator 
is also k-periodic. 

Proof. Consider the billiard inside the domain on the ellipsoid Mo bounded 
by the curve T A = M fl M A . We claim that the respective billiard ball 
map F\, acting on time-like geodesic segments, preserves the Joachimsthal 
integral J considered as a function on the space of geodesies. A similar fact 
for Euclidean ellipsoids is well known, see, e.g., [9| [22]. 

Consider a time-like geodesic segment 71 on the ellipsoid, reflecting in 
T A to another geodesic segment 72. According to Section [Tj, the straight 
lines, tangent to 71 and 72, are tangent to pseudo-confocal quadrics, say, 
M w and M^ 2 . We want to show that \ii = /i2- Indeed, the billiard system in 
the ambient space with reflection in the quadric Mx is integrable and M\ is 
orthogonal to M . Therefore the incoming and the outgoing rays are tangent 
to the same pseudo-confocal quadric, see Section (H and thus /ii = /x 2 . 
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According to Section [Q, F\ preserves the symplectic structure u> on the 
space of time-like geodesies; it also preserves the integral J. Therefore the 
1-form h(t)dt on the space of null geodesies is also invariant under the action 
of the map Fa on the space of null geodesies. 

In the limit A — > 0, the curve T\ becomes the tropic (T5]), and the billiard 
ball map F\, restricted to null geodesies, gets identified with the map T. 
Hence T preserves the 1-form h(t)dt. 

Finally, choose a cyclic coordinate s on the equator so that h(t)dt = ds. 
In this coordinate, the map T is a shift s i— > s + c. This map is /c-periodic if 
and only if kc G Z. This implies the statement of the theorem. □ 

Problem 5.2 It is interesting to find the relation on a, b, c, necessary and 
sufficient for the orbits of the map T to close up after k iterations and r turns 
around the equator. In the case of the Poncelet Porism, such conditions were 
found by Cayley, see [11] for a modern treatment. 

Remark 5.3 One also has Poncelet-style closure theorems for the geodesies 
in the polar caps and for space- or time-like geodesies in the equatorial belt 
reflecting from the tropics. Such results are similar to the ones known for 
Euclidean ellipsoid, see, e.g., [9j [22], and their proofs are similar to that of 
Theorem 15.11 but simpler: the T-invariant 1-form h(t)dt is obtained from 
a finite area form and a finite integral, not as a finite ratio of two infinite 
quantities. 

6 Curvature of the ellipsoid and a geodesi- 
cally equivalent Riemannian metric on it 

6.1 Curvature of the ellipsoid 

The behavior of the geodesies in a polar cap resembles that of the geodesies in 
the Poincare disc model of the hyperbolic plane; this observation is explained 
by the following proposition. 

Proposition 6.1 The Gauss curvature K of the ellipsoid is negative; it is 
given by the formula: 
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Proof. Similarly to the Euclidean case, one can use the normal Gauss map 
to compute the Gauss curvature of a surface, see [IE]. Consider the Northern 
polar cap. Normalize the normal vector as follows: 



so that (N, N) = —1. Thus the Gauss map sends the polar cap to the upper 
sheet of the hyperboloid of two sheets, and this map reverses the orientation. 
Hence the Gauss curvature is negative. 

Similarly, one considers the equatorial belt: the normal vector is then 
normalized to (N,N) = 1. 

It is straightforward but tedious to compute the Gauss curvature, and 
we do not dwell on this. The computation can be simplified by the use of 
Mathematica. □ 

6.2 Riemannian geodesically equivalent metric 

The next result follows from general constructions in [14j [19]; for complete- 
ness, we give a direct proof. Two metrics are called geodesically equivalent if 
they have the same non-parameterized geodesies. 

Proposition 6.2 The metric on the ellipsoid induced from the ambient 
Minkowski space, is geodesically equivalent to the Riemannian metric 



Proof. Let P = (u,v,w) denote a tangent vector at point Q = (x,y,z). 
Denote the diagonal matrix with the entries (1/a, 1 /b, 1/c) by A. To fix ideas, 
consider the equatorial belt. The Lagrangian for the metric (|T2|) is 




ds 2 




(12) 



L(P,Q) 



9 9 9 

U I I W 

a ~ l ~ b ~T c 



A(P) ■ P 

~1W~ 



x 2 I y 2 z 2 

? + F 7 1 



with 



a z cr 
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The Euler-Lagrange equation for a geodesic Q{t) with Q = P is 



L PP {Q)+L PQ {P)-L Q = \A{Q) 

where A(Q) is a Euclidean normal to the ellipsoid at point Q and A is a 
Lagrange multiplier; here Lpp and Lpq are the matrices of the second partial 
derivatives and Lq is the gradient vector. 
One easily computes: 

2 2 A( P) ■ P 

L " = /toy- 4 ' ipo = " m * v/w) ' L ° = " /4r v/(Q) ' 

and the Euler-Lagrange equation is rewritten as 

^-^p^ + w v/(g)=AA((?) - (13) 

To find the Lagrange multiplier, dot multiply equation (fT3l) by Q. One has: 

A(Q) ■ Q = 1, A(P) -Q = P- A(Q) = 0, V/(Q) • Q = 2/(Q), 
and hence 

A = A(g)-g + A(p)-p= ^ ( ^- g) =o. 

Thus (1131) implies that the acceleration Q lies in the plane spanned by the 
velocity vector P and the vector y4 _1 (V/(Q)) = 2N(Q), the Minkowski 
normal to the ellipsoid at point Q. It follows that Q(t) is a reparameterized 
geodesic of the restriction of the ambient metric on the ellipsoid. 
The case of the polar caps is similar. □ 



Remark 6.3 Note that the metric (fl2j) is given by the ratio of the two 
factors whose product is the Joachimsthal integral (jl]). 



7 Rigidity of ellipses as integrable Lorentz 
billiard curves 

In Section [331 we considered the limit of the ellipsoid (PQ) as b — *■ 0. In this 
limit, the ellipsoid becomes an ellipse in the Lorentz plane. This ellipse is 
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foliated by two parallel families of null lines, and the map T, defined at the 
beginning of Section El takes oriented lines in one null direction to oriented 
null lines in the other one. 

More generally, consider an oval (closed smooth strictly convex curve) 7 
in the Lorentz plane and the billiard system inside it. The billiard reflection 
takes one null direction to the other. Thus, restricted to the null lines, the 
billiard ball map is identified with the circle map T : 7 — ► 7 described 
as follows. The intersections with the lines in the first direction define an 
involution on 7, and likewise for the second direction. The composition of 
these involutions is the map T. 



The circle map T is depicted in figure [TOl it was previously considered 
in different contexts: in relation to Hilbert's 13th problem [3]; the Sobolev 
equation, approximately describing fluid oscillations in a fast rotating tank, 
[T8] ; and the theory of Lorentz surfaces, where the rotation number of T 
provides a continuous invariant of the conformal class of a Lorentz disc [23J . 

The map T : 7 — ► 7 depends only on the choice of two (null) directions, 
say, u and v; thus we denote it by T( U)V y If 7 is an ellipse then T( u ,v) is 
conjugated to a rotation. This is a limit case of Theorem 15.11 but it is also 
easily proved directly: an afline transformation takes 7 to a circle, and then 
T( u ,v) becomes the rotation through the angle twice that between the two 
respective directions. 

Conjecture 7.1 Let 7 be a plane oval such that, for every pair of directions 
u and v, the map T( Ujl) ) is conjugated to a circle rotation. Then 7 is an ellipse. 

We prove a weaker result. 




Figure 10: A map of an oval 
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Theorem 7.2 Let j(t) be a parameterized plane oval with the property that, 
for every pair of directions u and v, the map T is a translation in the variable 
t, i.e., there is a constant c(u,v) such that T( U)V )(y(t)) = j(t + c(u,v)). Then 
7 is an ellipse. 

Proof. Without loss of generality, assume that 7 is parameterized by t G 
R/2ttZ. 

Choose a direction u. There is a unique maximal chord of 7 in direction 
it, say, AB; such a chord is called an affine diameter. Since AB is maximal, 
the tangent lines to 7 at A and B are parallel; let v be their direction. It is 
convenient to apply an affine transformation that makes u horizontal and v 
vertical. 

Let C and D be the points at which the tangent line to 7 are horizontal. 
We claim that the chord CD is vertical. 

Indeed, consider the map T( u>v y Then T 2 (A) = A. Since T( u>v ) is conju- 
gated to a rotation, every point of 7 is 2-periodic. Let CD' be the vertical 
chord of 7. Since T?JC) = C, the direction of 7 at D' is horizontal, and 
hence D' = D, see figure [TT1 



c 




D' D 



Figure 11: Proving that D = D' 

Thus, for every direction u, there exists a conjugate direction v: there are 
affine diameters AB and CD of the curve 7, having directions u and v, such 
that the tangent lines at A and B are parallel to CD, and the tangent lines 
at C and D are parallel to AB (this makes 7 a Blaschke P-curve [7]). 

Furthermore, if (u, v) is a pair of conjugate directions then every point of 
7 is a vertex of an inscribed parallelogram with the sides parallel to u and v, 
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and whose diagonals correspond to a 2-periodic orbit of T( UjV y The opposite 
vertices of these parallelograms are 2-periodic points j(t) and j(t + 7r). 



















Figure 12: Inscribed rectangles 

Assume that 7 is centrally symmetric with respect to the origin O. Then 
each rectangle is centered at O - otherwise 7 is not convex, see figure [T2J 
It follows that the mid-points of the horizontal sides of the rectangles lie on 
the affine diameter CD, and the mid-point of the vertical sides lie on AB. 
This is true for every pair of conjugate directions, hence 7 admits affine line 
symmetry for every direction. This is a characteristic property of ellipses [6] , 
and thus 7 is an ellipse. 

It remains to prove that 7 is centrally symmetric. Consider two directions, 
u and v, making an infinitesimal angle e. Consider figure [TBI in which AB is an 
affine diameter in the direction u, and the lines BC and AD have the direction 
v. One has: T (UjV) (A) = C, T (u>v) (B) = D. Let A = 7(f), C = ^{t + 5). Then 
B = j(t + 7r), D = j(t + ii + 5). Since the opposite sides of the infinitesimal 
quadrilateral ACBD are parallel, it is a parallelogram and AC = BD. It 
follows that j'(t) = — 7 / (t + 7r), and hence j(t) +j(t + vr) is a constant vector. 
Thus 7 is centrally symmetric. □ 

Acknowledgments. We are grateful to V. Chernov, D. Fuchs and P. Lee 
for stimulating discussions. The second and the third authors were partially 
supported by NSERC and NSF grants, respectively. 
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B D 

Figure 13: Proving central symmetry of 7 
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